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ABSTRACT. The expression of the solutions of a special boundary value prob- 
lem determined by an ordinary linear n-dimensional differential equation is 
given when the deviations are mixed and they have fixed points. The key 
assumption is that the deviations are odd functions. 


1. INTRODUCTION 


In this note we deal with the linear n-dimensional functional differential equation 
of mixed type of the form 


(1.1) DEaps: i(t) = A(t)x(t) + B(t)2(t — r(t) + f(t), tER, 


where 7 is a real valued (nontrivial) odd continuous function. Thus, if in a set J, 
the function 7 is positive we have delay argument. Then in the symmetric set —I 
it is negative and so we have advanced argument. The motivation of this paper 
comes from [2] due to L. Berezansky, E. Braverman and S. Pinelas, where, in the 
first one, a differential equation of the form 


a(t) + $ dg (t)a(hg(t)) =0, t> to 
k=l 


is investigated. (For a partial case see [3].) Here the functions bẹ and t > t — 
h(t) are Lebesgue measurable and essentially bounded on an interval of the form 
[to, +oo) and satisfy —o, < t — hk(t) < Tk, with on, Tk > 0. (See, § 1.4 in [2].) In 
that paper the authors are interested in the existence and uniqueness of solutions, 
as well as, in the asymptotic properties of the solutions. Our interest is focused on 
their remark that for the mixed equation 


(1.2) é(t) + b(t)x(t — sint) = f(t) 


the authors are not aware of any existence and stability results. (Relation (1.1) 
in [2].) Also, at the end of the paper the authors among others in a discussion 
of the results (no. (2)), pose the following open problem: Is it possible to find 
conditions such that there exists a unique solution (not necessarily bounded) of an 
initial value problem for some types of equations with mixed arguments. In this 
paper we present such a type of differential equations and give answers to these 
questions. To be more precise we formulate a initial value problem and give the 
expression of the (unique) solution in the form of a series with terms containing 
iterated integrals. To obtain our results we use the simple (but very important) 
fact that the roots of the argument 7 are fixed points of the function t > t — T(t). 
Then the results of [7] play a very important role. 
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To shorten this work, we do not intent to go back to refer to the history of the 
mixed type differential equations, but, instead, we can say that more facts about 
differential equations with mixed arguments can be found in a great number of 
works in the literature, starting from the work due to Kato and McLeon [8]. See, 
also [1, 2, 3, 4, 5, 9] and the references therein. 


2. THE MAIN RESULTS 


We start with some facts borrowed from [7]. Let J := (a,b) be an interval of the 
real line R. Consider the linear n-dimensional delay differential equation 


i) = A(t)a(t) + Bit)x(g(t)) + FE), 


where 

(C1) g: I > I is a continuous function such that g(t) < t, for all t € I and for 
some w with g(w) = w it holds g(t) > w for t € (w,b). 

(C2) A,B are n x n-matrix valued continuous functions defined on J and inte- 
grable on (w, b). 

(C3) f : I> R” is a function continuous and integrable on (w, b). 

We furnish the space of n x n matrices with a norm || - || induced from the 
euclidean norm of the vector space R”. 

From the standard theory of linear functional differential equations (see, e.g., 
[6], p. 11) it is well known, that under these conditions, for any o € I and any 
continuous function y : [9(07),0] —> R”, there exists exactly one solution x(t) := 
x(t; (o, p)) such that x(t) = y(t), for all t € [9(0), o]. Here g(o) = inft>o g(t). Let 
w be a fixed point of g on I; then [g(w), w] = [w] and therefore any initial function 
y for the solution x(w, vy) is actually a point of R”. Also, if g does not have a fixed 
point w in J, as, for example, we have in case of constant delay, namely, g(t) = t-r, 
but it happens lim;_,,,4 g(t) = w E€ RU {—oo}, then we say that w is a generalised 
fixed point of g. In this case the solution x, whose limit lim;_,.,4 x(t) =: (w) = € 
exists, it is unique and it is given by 


x(t) = Ey p(t,w)€+ z(t), tE (w,d), 


where £4 p(t,w) is the evolution matrix of the homogeneous equation and z is the 
solution of the non-homogeneous equation such that z(w) = 0. Here the evolution 
E4,g is defined by the formula 


+00 
(2.1) Ea.a(t,w):= X ([Z2(A, ByInxnl(t), if t € (w,d), 
m=0 


where Inx» is the n x n identity matrix and 


R(t), if m=0 
[Z2 (A, B)R\(t) = 4 Ji [A(WDIZZ (A, B)R](u)+ 
+B(u)[Z®-1(A, B)R\(g(u))]du, if m> 1, 


for all continuous n x n-matrix valued functions R defined at least on the interval 
(w,b). Notice that the evolution satisfies the equation 


S Ea p(t,u) = Al Ea.pltw) + BOBA Bll) u), tE (wb). 
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Also, the solution z of the non-homogeneous equation with initial value z(w) = 0, 
is given by 
+00 
2(t) = $ [ZIA B)F]®), 
m=0 
where F(t) := fi f(u)du, t € (w,b). 

In the present case the basic assumption on the argument 7 is the following: 

(C4) The function T is continuous and odd and such that in any compact interval 
there is a finite number of points where it changes its sign. 

By assumption (C4), the roots of 7 are distinct. Let rj, j € {0,1,---,N}, 
(N € N), be the positive moments where the lag 7 changes sign. Hence, due to 
oddity of 7, the points —r,;, j € {1,2,--- , N}, are the negative zeros of 7 the points 
—rj,j € {0,1,--- , N} are also zeros of the lag. We shall denote by J,, the interval 
[rj,7j41], for all j € Z. These items are such that, if in some interval J,,, it holds 
T(t) > 0, then in the next interval [,,,,, it holds T(t) < 0. Moreover, if it holds 
T(t) > 0 on the interval J,,, then T(—t) < 0 on the symmetric interval J_,,,, and 
inversely. In other words, in this case, the equation is a delay differential equation 
for the interval Ir and an advanced differential equation in the interval [_ 

We, further, assume the following: 

(C5) The function g(t) := t — T(t) maps any I,, to itself, for all j € Z. 

The prototype of such a function 7 is the trigonometric function sin used in (1.2). 


Tj+1' 


T5 r4 T T3 T2 T rı 0 + rı T2 T T3 r4 T T5 T6 


In the intervals where 7 gets positive values we have delay differential equation, while in the 
intervals with negative values of r we have differential equation of advanced type. 


To proceed we need the following lemma: 


Lemma 2.1. Assume that conditions (C1, C2, C3, C4, C5) are satisfied and let 
I := |a, 8] be an interval of the form I,,, j € Z. Then, for any € € R”, there exists 
a unique solution of (1.1), defined on I, with the initial value x(a) = €. Moreover, 
for any Å ER, there exists a unique solution of (1.1), t € I, with the terminal value 


(8) =¢. 


Proof. The proof goes as in Theorem 2.1 in [7], provided that we choose tg such 
that 


max{to,g(to)} 
(2.2) 1 (A)| + I|B(u)l])du < 1. 


The second result is proved by using the first part and the fact that by the 
transformation t + —t =: s, a terminal value problem concerning a differential 
equation DE, p,¢ becomes a initial value problem of the form 


(2.3) DE; gf: 9(s) = A(s)y(s) + B(s)y(9(s))) + f(s), s € [-B,—al], 
with y(— 8) = ¢, where 
(2.4) A(s) =—A(-s), B(s)=—B(-s), f(s) = —f(-8). 
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From the arguments above, it seems that an initial value problem formed by the 
differential equation (1.1), associated with the initial value z(w) = €, behaves like 
the analogous one for linear ordinary differential equations. And this is true because 
the initial value is a vector and not a function defined on the ” past”. Moreover, 
in the ODE case the corresponding evolution is produced by the transition matrix 
and it is well known that it is a nonsingular matrix. In the present functional case, 
this is not true, as it is proved in the scalar case in [7]. The question is under what 
conditions the evolution is a regular matrix. In the real-scalar case easy conditions 
can guarantee this fact as follows: 

Consider an interval [a, 2] of the form J, and the homogeneous part of equation 
(2.1). Also, assume that the functions A, B are both real nonnegative, odd and 
continuous, or they satisfy the condition! 


Tj+1 
(2.5) J (A+ [Baus m2, j= 0,15. 
Then, in case T(t) > 0, from (2.1) we get 

Ea p(t, a) = I; te la, 8], 
in the first case, and 

+00 1 


Ealt) 21- $ C (Awl + BOD)” 


m=1 ‘ 


t 
= 2— exp | (|A(u)| + [B(u)|)du > 0, t€ [aA] 
in the second case. Also, in case T(t) < 0, from (2.3) we have 
E4, B(s, a) 2 1, SE [-8, —a], 


in the first case, and 
+oo 


1 
Ea p(s,a) >1- 5 mi 


=2- exp f (\AQW)|+|Blu))du> 0, s€ [2,0], 


in the second case. Notice that it holds Â = A and Ê = B. Thus, in all cases, the 
evolution is not zero. 
For the n-dimensional case we give the following result: 


Cf (Aol + Beda)" 


Proposition 2.2. Let |a, 8] represent an interval I,,, such that T(t) > 0, for 
t € (a, 8). Assume, also, that the matriz-valued functions A,B are odd functions. 
Then the evolution E4 g satisfies the identity 


Ea,p(—t, —B)Ea,3(8,a) = Ea, g(t a), tE [a, £] 
and therefore, we have 
Ea p(-a, —6)E4,B (B, a) = Inxn- 


The latter implies that the matrix E4 B(6, œ) is regular and its inverse is the matriz 
E4, B(—a, —B). 


Notice that, in this case, relation (2.5) implies (2.2). 
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Proof. First notice that Â = A and Ê = B. Also, the points —8,—a are fixed 
points of g(t) = t — T(t). Let x be a (nontrivial) solution of the homogeneous 
part of equation DE 4,p,7 on the interval [a, 6]. Thus the unique solution y of the 
homogeneous part of equation (2.3), with initial value y(—8) = x(3), can be written 
as 
y(s) = Ea,a(s,—B)y(—B) = Ea,a(-t,-8)2(8), tE [a, 8]. 

Define the function w by w(t) := y(—t), t € [a, 6]. Since the point —£ is fixed for 
g and Ez is the evolution of the equation, we have 


.B(—t, —6)a(8) + B(t)Ea,p(9(-t), —8)a(8) = 
=t) + Bly g(t)) = A(t)w(t) + Bit)w(g(t)). 
) 


But we have «(8 a w(8) and so, by the uniqueness of solutions (Lemma 2.1), i 
follows that, for all t € [a, 8], it holds 


0 = z(t) — w(t) = z(t) — y(-t) = Ea p(t, a)x(a) — Ea,p(-t, —8)x(8) 
= [Ea,p(t, a) — E4,B(—t, —8)E4,3(8, @)] x(q). 


Since x(a) could be any vector in R”, it follows that the null space of the matrix 
in the brackets is the whole space R” and so we must have 


Es p(t, a) = E4 p(-t, —B)E4 p(B, a) = Onxn: 
which proves the result. 


“ | 
= 
= 
x 
€ 


Our main result of this work is given in the following theorem: 


Theorem 2.3. Consider the mized type differential equation (1.1) and it is assumed 
that assumptions (C1, C2, C3, C4, C5) are satisfied. Assume, also, that the matrices 
E4.p(—reai-1, —rat) are nonsingular, for all l = 1,2,3,---. Then, for any € € R”, 
there exists a unique solution x(t), t E R such that x(0) = £. 


Proof. Our schedule is to apply the idea of [7] to each interval of the form Ir, 
and J_,,,, and obtain successively the expression of the solutions on each such 
interval. The problem is what happens when in an interval of the form (a,b) it 
holds t — T(t) > t, namely the differential equation is advanced. In order to face it 
we use the symmetric interval (—b, —a), where we have delay argument. 


a a aT 


Ts T4 T T3 T2 T rı 0 + rı T2 T T3 T4 F T5 T6 


Our method relies on the following procedure: If in an interval we have delay 
argument we apply the method described in [7] and, if in an interval the argument 
is advanced, then by the transformation t + —t we use its symmetric one, where 
we have delay and then we come back to the original interval. This justifies our 
assumption on the oddity of 7. Finally, expressing the solutions on each interval 
piece-wise we ”join” them by using the initial values at the fixed points of the 
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argument, namely at the terms of the sequence (rj). To be more precise, we solve 
the problem step-by-step as follows: 

Step 1. Solution on Ip: First we assume that it holds T(t) > 0 on the interval 
Io. (If T(t) < 0 on this interval we start from the interval J_,, and follow analogous 
steps.) Then we have the delay differential equation DE4,g,f, namely the equation 


(2.6) a(t) = A(t)a(t) + B(t)x(t — r(t)) + f(t), tE Io. 


Since 0 is a fixed point of the argument g(t) := t— T(t), according to [7], the solution 
of (2.6) with initial value € at 0 is uniquely determined by its initial value € at 0 
and it is given by 


+00 
(2.7) x(t) = Ea, B(t,0)E + XO [Z0 (A, B)Fol(t), t€ Lo, 
m=0 
where Fo(t =f F 
R(t), if m=0 


[25°(A, B)R]() := 4 fo [AQ [29° (A, B)R\(u) + 
a [Ze"(A, B)R\(g(u))]du, if m>1 
and E4 p(t,0) is the function given by (2.1) for w = 0. 

Step 2. Solution on I_,,: On this interval it holds T(t) < 0 and so g(t) := 
t— T(t) > t, ie. equation (2.6) is advanced. Observe that for t € I_,,, we have 
s = —t € Ip. So, making the transformation y(s) := x(t) and (2.4), we get the 
equation 


(2.8) DE; gz: UCs) = A(s)y(s) + B(s)y(g(s))) + f(s), 8 € To, 


which is a delay differential equation, due to the fact that it holds g(s) = s—rT(s) < 
s, for all s € Io. Then, as above, we can solve it and obtain its solution 


+00 
y(s) R Ex g(s,0)£ + 3 [Z5 (A, B) Fo](s), sE Io, 


m=0 


where Fy(t =f f u)du. Notice that y(0) = €. Now, let w(t) := y(—t), t € I-n. 
Since Ey. is athe a of equation DE, p f we can easily get that 


b(t) = —y(-t) = —A(—t)y(—t) - B(-t)y(9(-1))) — f(t) = 
= A(t)w(t) + B)w(g(t))) + FO), 


where, notice that y(0) = 7(0). Therefore x = w and so, the unique solution x of 
the original problem on the interval J_,, is given by 


(2.9) a(t) = —t, 0)E + 3 Z? (Â, B)Ê](—t), t€ In. 
m=0 
Step 3. Solution on I_,,: On the interval I—,, equation (2.6) is a delay differ- 
ential equation and —rg is a fixed point of the delay. Hence, its solution with initial 
value z(—r2) is given by 
+00 
(2.10) a(t) = Ea.p(t,—r2)2(—re) + X [Z™,,(A, B)F-r](t), tE Longs 


m=0 
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where F_,.(t = f ra f (u)du and 
(2.11) 
R(t), if m=0 
[27,.(A, B)R\(t) = 4 JÉ, [A(w)[Z-r2(A, B)R] (u)du, 
+B(u)[Z™,"(A, B)R](g(u))] du, if m> 1. 


—T2 


Our purpose is to express the value «(—r2) in terms of the initial value €. From 
(2.10) we have 


+00 
w(—r1) = Ea,p(—?,—?2)@(—re) + y AEB N 
m=0 
and from (2.9) 
a(—1r1) = E; p(r1,0)§ + z [Z0 (Â, B) Fo](r1). 


These two equations give 
-1 
x(—r2) = (E4,B(-r1, —r2)) one 


1 = ZE (Â, Ê)Fo](r1) — (2%, (A, B)F-ra](=r1))]. 
Therefore on the interval 7 -r> the solution x has the expression 


z(t) = Ea,p(t,—r2)(Ba,a(—r1, —r2)) [Bap (rn, O)E+ 


+ pak ZE (Â, B)Fol(ra) — (2, (A, B)F-ral(=r1))| + 
(2.12) + Suen m (A, B)E rll), tE Lrg 
Step 4. Solution on I,,: On this interval the argument g(t) = t — T(t) is 
of advanced type and, for t € I}, we have s = —t € I_,,. So, making the 
transformation (2.4) and by setting y(s) = x(—s) the equation becomes 
(2.13) g(s) = A(t)y(s) + B(s)y(s — r(s)) + f(s), 8 € Lone, 


which is a delay differential equation, due to the fact that it holds g(s) = s—r(s) < 
s, for all s € I_,,. But from (2.12), we know that such an equation has the solution 
expressed as 


y(s) = (Ea,B(=r1, —r2)) "EA, B(s, —r2) |e a B(r1,0)£+ 


+ > (Z (Å, Ê)Fol(ri) — 12%, (4, B)F-ra](—r1))| + 
m=0 
+00 
(2.14) + X [2™.,(A, B)F-ra](8), t € Ir. 


m=0 
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where F_,.(s og u)du and Z™.,(A, B) is given by (2.11). Therefore on the 


interval Z., É e = ) = x(-; (0, £)) has the expression 


E een oe eer 


+ (IZE (Â, B) Fol(r1) = [Z (A, B)P_ra] (ri) | + 


(2.15) + Salam ‘(ABP 3 ). t€ In. 


Step 5. Solution on L: a proceed to the interval J,,. In this interval the 
argument is a delay, so the solution with initial value x(r2) at r2 is given by 


(2.16) a(t) = Ea.p(t,r2)x(re + Siena. aR mlt), tE Ln, 
where F,,(t =f f(u)du, and 
R(t), if m=0 


[272 (4, B)R]() = 4 Jo [AWZ (A, B)R](u)+ 
uZ (A, B)R](g(u))]du, if m > 1. 
The value z(r2) can be found from (2.15) for t = rg. Hence the solution on [,.. is 


given by 
x(t) = Ea Blt, r2) — —r)) |E4,a(r O)E+ 


+ x Z (Â, B)Fol(rs) — [Z2 (4, B)F-ral(—r1)) |+ 
(2.17) 
+00 
+ yen ot As a — ra)) + XC [ZR (A, B)Frol(t), tE Ira, 
m=0 


because E4,B(—r2, a = Tagi 

So far we have solved equation (2.6) on the entire interval [—r2, r3], where we gave 
the solution on the subintervals Jo, [_,,, [-r., Ir,, Ir, and we expressed the unique 
solution respectively by the types (2.7), (2.9), (2.12), (2.15) and (2.17). Now, we 
can continue in the same process on the intervals I_,,, I-,, Ira; Ira; Drs, Lrg; 
Irs, Ire, etc. and inductively on the intervals I-ra; Loroj42)1 Irez Trea? 
for all integers j and obtain the expression of the solution «x(-; (0, €)) in any interval 
of the real line. The proof is complete. 


3. AN APPLICATION 


As an application we borrow an example from [7], considering, for simplicity, 
only the homogeneous part of the equation. Consider the scalar linear differential 
equation 

ilt) = at’a(t — pt) +t?, tER, 
where a,b > 0, d > 0 and p € (0,1). Here the function 7(t) := pt, t € R is odd with 
only one zero t = 0 and it satisfies the condition (C5) on the symmetric intervals 
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(—oo, 0) and (0, +00). Its only zero is 0. Let c := 1 — u (€ (0,1)). Applying the 
previous arguments presented in Theorem 2.3, step 1 and step 2, we obtain? that 
the solution x with initial value z(0) = € is given by 


(3.1) 

+00 a kg k(k— 1)(b+1) k(b+1) 
k=0 REDE | ft 

¿9+1 Lyo akck(d+1)+4k(k-1)(+1) k(b+1) t>0 

d+1 k=1 (otd2)(PtdE3) (ROGET) eZ 

x(t) = a 
(t) MARS ak eR RR 1)(b a tree, 

k=0 kI(b+1)* 

E pati (= a)*e k(d+1)+$ k(k—1)(b+1) 


aa |1 eis 1 Grice aon TEFA FTT ( geo, t<0. 
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